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Abstract 



Studies of quarkonium spectral functions at finite temperature, based on an approach combining QCD sum rules and the maximum 
entropy method are briefly reviewed. QCD sum rules for heavy quarkonia incorporate finite temperature effects in form of changing 
values of gluonic condensates that appear in the operator product expansion. These changes depend on the energy density and 
pressure at finite temperature, which we extract from quenched lattice QCD calculations. The maximum entropy method then 
allows us to obtain the most probable spectral function from the sum rules, without having to introduce any specific assumption 
about its functional form. Our findings suggest that the charmonium ground states of both S-wave and P-wave channels dissolve 
into the continuum already at temperatures around or slightly above the critical temperature T c , while the bottomonium states are 
less influenced by temperature effects, surviving up to about 2.5 T c or higher for S-wave and up to about 2.0 T c for P-wave states. 
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^ . 1. Introduction 

The study of quarkonia in hot matter has, since the early suggestions (jl that were made more than 25 years ago, 
evolved into a field with diverse activities in both experiment and theory. Especially, through the heavy-ion collision 
experiments at RHIC and LHC, a huge amount of experimental data on quarkonium production in various reactions 
is now available, which can be compared with theoretical expectations. This task, however, has turned out to be very 
complex, as a large number of competing effects have to be taken into account to describe the experimental results (for 
reviews see JUSl). The most basic inputs of these calculations are the quarkonia spectral functions, which include 
all the physically relevant information of the quarkonium states as well as their behavior at finite temperature. It is 
therefore desirable to obtain these spectral functions from calculations based on QCD. Progress in this direction has 
been acheived mostly due to the advancement of lattice QCD and the use of the maximum entropy method (MEM) for 
the extraction of the spectral function from the euclidean time correlator J3l (see also Jit and the references therein 
for the latest results). 

Another approach for capturing information on quarkonia spectral functions at both zero and finite temperature is 
provided by QCD sum rules. This method exploits the analytic properties of the two-point function of operators to 
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derive certain integrals over the hadronic spectral functions (the "sum rules"), which, via the operator product expan- 
sion (OPE), can be related to a combination of perturbatively calculable quantities and non-perturbative condensates, 
containing information on the QCD vacuum. In the case of the quarkonia channels considered here, these are gluonic 
condensates, the most important one being the gluon condensate of mass dimension 4 

Recently, it has become possible to make use of the MEM technique to analyze QCD sum rules |7|], which allows 
to extract the most probable form of the spectral function from the OPE without having to resort to some specific 
functional form. This approach has since been applied to both charmonium 0] and bottomonium |9] channels and to 
discuss these results will be main goal of this article. 

The proceedings are organized as follows. After a brief description of our formalism and analysis methods, we 
recapitulate the results obtained so far in J§,0]. We will also show some novel results on the thermal behavior P- 
wave charmonia (xm, Xc\), which were not included in |@]. Finally, we summarize our results and give an outlook 
by discussing open issues and possible future directions of further improvements of the sum rule approach presented 
here. 

2. Formalism 

QCD sum rules at finite temperature make use of the analyticity of the two-point correlator of some general local 
operator f(x): 



Here, f(x) stands for hy^hix), hysh(x), hh(x) and (q^qy/q 2 - g^v)hy5y v h(x) in the vector, pseudoscalar, scalar and 
axial-vector channel, respectively, while h(x) represents either a charm- or bottom-quark field. The definition of the 
expectation value (0)t is (0)t = Tr(e~ H t T 0)/Tr(e~ H/T ). Through a dispersion relation, one can relate the correlator 
calculated in the deep-Euclidean region (-q 2 — > oo) to a specific integral over the hadronic spectral function p J (s). 
After the application of the Borel transform, one arrives at the following expression: 



In this equation, the left-hand side can be calculated analytically using the OPE. To get the spectral function p J (s) 
one therefore somehow has to invert the above integral. This is, however, an ill-posed problem and can not be solved 
rigorously, because the left-hand side is only known as an asymptotic expansion, with coefficients determined with 
limited precision. With the help of Bayes' theorem, it is nevertheless possible to obtain the most probable form of 
p J (s), given Eq. (0 and additional information on the spectral function such as positivity and asymptotic values at high 
and low energy. This is the essence of the MEM approach discussed in this article. For the concrete implementation 
of MEM for analyzing QCD sum rules, see |7|]. 

For the sum rules of quarkonia, all the finite temperature effects can be included into the temperature dependent 
condensates owing to the large separation scale in the OPE lioll . This speration is valid as long as the temperature 
does not reach values of the order of the quark mass. The temperature dependencies of the condensates have been 
obtained in [5], where the relation between the gluon condensate and the energy momentum tensor is exploited to give 
the value of the gluon condensate as a function of the ener gy d ensity and pressure. These thermodynamic quantities 
are then extracted from quenched lattice QCD calculations 111 ill 1211 . In this way, we are able to calculate the left-hand 
side of the sum rules as a function of temperature. As a last step, we then use MEM to retrieve the corresponding 
spectral functions from Eq. (0. 



Let us first discuss the results of the charmonium channels. They are shown in Fig. Q] The S-wave channels (left 
plots) have already been presented in J3l, while the P-wave channels are shown here for the first time. Concentrating 
firstly on the spectral functions at zero temperature, it is seen that clear peaks are generated, which represent the 
lowest state of each channel. The positions of these peaks reproduce the experimental values with a precision of about 
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Figure 1. Charmonium spectral functions at zero and finite temperature in the pseudoscalar (top left), sclalar (top right), vector (bottom left) and 
axialvector channel (bottom right). The left figures are adapted from @|. 



At finite temperature, we observe that the lowest peaks of all channels vanish slighly above the critical temperature 
T c . The origin of this melting effect is a sudden change of the gluonic condensates around T c , which can be related to 
the deconfinement transition of the gluonic matter. 

Next, let us look at the results for bottomonium. These are given in Fig. [2] Here, as for charmonium, clear peaks 
are seen for each channel at zero temperature. These peaks are found at 150-500 MeV above the experimental values 
of the respective ground states. This discrepancy is caused by the excited states (for instance T(25) and Y(3S) in the 
vector channel), which cannot be resolved by the MEM analysis, and pull the lowest peaks to higher energies than the 
actual ground state. This is in contrast to the charmonium case, in which the excited states give only a relatively small 
contribution to the lowest peak. For a more detailed discussion of this issue, see |@]. 

Turning now to the finite temperature curves, it is noted that the bottomonium states are modified much slower 
than their charmonium counterparts, which is in agreement with phenomenological expectations. Concretely, both 
S-wave spectral functions still exhibit a clear peak at T — 2.0 T c which starts to dissolve at about 2.5 T c . The P-wave 
states are on the other hand modified somewhat faster and disappear already at temperatures around 2.0 T c . The reason 
for the robustness of the bottomonium states can be traced back to the fact that the gluon condensate terms in the OPE 
are proportional to 1 /mf , m/, being the quark mass. These are the driving terms of the quarkonium melting, and are 
therefore relatively suppressed for the bottomonium sum rules. This is why for heavier quark masses one needs to go 
to higher temperatures to observe a significant effect. 



4. Conclusions and Outlook 

We have analyzed quarkonia ground states at finite temperature in the vector-, pseudoscalar-, scalar- and axialvector- 
channels. By combining the techniques of QCD sum rules and MEM, we have extracted the spectral functions from 
the OPE of the correlators calculated in the deep Euclidean region. As a result, it is found that the charmonium ground 
states of both S-wave and P-wave channels dissolve into the continuum already at temperatures around T c , while the 
bottomonium states are more stable, surviving up to about 2.5 T c or higher for S-wave and about 2.0 T c for P-wave 
states. 

On a qualitative level, the results shown in the previous section appear to be quite reasonable. We, however, have 
assumed that the OPE converges sufficiently fast and thus the higher order terms will consist of only small corrections 
by which the general picture is not altered. The only way of checking the validity of this assumption, is to take into 
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Figure 2. Bottomonium spectral functions at zero and finite temperature in the pseudoscalar (top left), sclalar (top right), vector (bottom left) and 
axialvector channel (bottom right). All shown figures are adapted from @|. 



account more higher order terms and to evaluate their effect on the spectral functions. The OPE generally contains 
two sorts of terms: perturbative and non-perturbative. The perturbative terms include the higher order a s corrections 
to the Wilson coefficients appearing in the OPE. Among them, especially the second order correction to the identity 
operator Wilson coefficient is potentally large 11311 . and a detailed calculation of its contribution to the sum rules of 
the vector channel is presently ongoing [14]. Furthermore, the non-perturbative corrections contain various gluonic 
condensates of higher dimension, for whose one ideally needs to know the value at zero and finite temperature. Such 
information can only be accessed via a lattice QCD calculation, which is planned to be carried out in the future. 
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